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1. Introduction 

Today most lattice QCD simulations are performed at unphysical light dynamical quark masses. 
Chiral perturbation theory (ZPT) [[I], ^ has proved to be a very important tool for such simulations. 
By using ZPT, one can extrapolate physical quantities to physical light quark masses and get in- 
formation on low energy constants (LECs) of the chiral theory. Although three-flavor ZPT has 
been successfully used for simulations with 2+1 dynamical quarks, we are still interested in the 
two-flavor ZPT for the following reasons: 1) Usually the simulated light quark masses are much 
smaller than the simulated strange quark mass. We expect the SU(2) expansion to serve as a better 
approximation and to converge faster than the SU(3) one. 2) Fits to SU(2) ZPT can give us di- 
rection information about LECs in the two-flavor theory. 3) By comparing results for SU(2) and 
SU(3) fits, one can study the systematic errors from truncations of different versions of XYI. 

In this work, we study the SU(2) ZPT for staggered fermions in the partially-quenched case, 
and obtain relations between SU(2) and SU(3) LECs by comparing formulae for the pion mass and 
decay constant from SU(2) and SU(3) ZPT. Then, we perform a systematic NNLO SU(2) chiral 
analysis for recent MILC data in the light pseudoscalar sector. Results for the pion decay constant, 
SU(2) LECs and chiral condensate in the two-flavor chiral limit are presented. 



2. Rooted SU(2) staggered chiral perturbation theory 

For lattice simulations based on the staggered fermion formalism, the correct effective theory 
is rooted staggered XVY (rSZPT) g, Q], in which taste- violating effects at finite lattice 

spacings are incorporated systematically. Physical quantities expressed in rSZPT become joint 
expansions in both m q and a 2 , where a is the lattice spacing. The three-flavor rSZPT has been 
well established [Q] and successfully applied in analyzing the lattice data. Here we concentrate 
on the two-flavor case. Instead of the three-flavor chiral limit m u = = m„ = 0, we perform the 
expansion around the two-flavor chiral limit m u = m d = 0,m s = m^ hys ', where m p s hys is the physical 
strange quark mass. 

The SU(2) rSZPT can be constructed by following the same procedure used for SU(3) rSZPT [Q]: 
First, one writes down the Symanzik effective theory (SET) for staggered fermions. Second, one 
maps the terms in the SET to operators in the chiral Lagrangian by using a spurion analysis. The 
power counting rule depends on the specific version of staggered fermions being used. For asqtad 
staggered fermions, we use @,|] a 2 8 ~ 2Bm ~ p 2 where 8 is a typical taste-splitting term. 

At leading order (&(a ,p ,m q )), the chiral Lagrangian for SU(2) SZPT is 

j2f(4) = ^Tr(D M ED M Et) _ ^Trfctf + X L) 

+ ^f-(U} l + ... + U? n > + D) 1 + . . . + D^) 2 + a 2 f , 
X = 2H(2)Diag( m x I,... ,m x I ,m y I, . . . , m y I, m u I, . , m u I , m d I, . . . , m d I ) , (2.1) 

ri T n ' r "r "r 

where E = exp (/<!>//) and ~f is the LO taste-violating potential. Their definitions can be found 
in Ref. [Q]. In Eq. ( |2.1[ ), the replica method is used explicitly: we take n' r copies of each valence 
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quark and n r copies of each sea quark. At the end of the calculations, we set n' r = to account for 
partial-quenching, and n r = 1/4 for taking the fourth-root of the fermion determinant. 

At NLO, the SZPT chiral Lagrangian contains two parts: the continuum terms at order £?(p 4 ,p 2 
m q ,m 2 ) and taste-violating terms at order &(p 2 a 2 ,m q a 2 ,a 4 ). In the partially-quenched case, the 
NLO continuum SU(2) chiral Lagrangian reads: 

■S& = "f [Tr^Z 1 ^!)] 2 - |Tr(D M I t D v £)Tr(D A1 I t D v £) 
+ ^(Tr(D A1 £ t D M £D v £ t D v £) - ^[Tr(D M EtD M E)] 2 
+ pi (Tr(D A1 I t D v £D /J £ t D v I) + 2Tr(D /J I t D M ID v I t D v I) 

- ^[Tr(D^D^)] 2 -Tr(D M £ t D v I)Tr(D /i £ t D v r 

- ^[TrCfEt + E^t)] 2 + |Tr(D M lt DM£)Tr(z£ t + £z t) 

+ & (zTrp? z^Z + *zW) ~ MX^ + ~ Tr(z^ - 
+ |[Tr(^_ £z t )] 2 

if 

- l° 5 Tr(tfF R ^LF L ^) - -j-Tr(F LtiV D^D v L + F Rflv D^D v ^), (2.2) 



where are bare LECs in ordinary SU(2) ZPT, and p®-p® are four new (bare) LECs in SU(2) 
PQZPT. Operators associated with p® are unphysical operators, in the sense that physical matrix 
elements of these operators vanish in the limit where the valence quark masses are set equal to sea 
quark masses. The specific form of the NLO taste-violating terms are not relevant to this work, so 



we do not list them here. For more details, please see Ref. [10] 



With the LO and NLO SU(2) SZPT Lagrangian, one can calculate the pion mass and pion 
decay constant in the partially-quenched case. Throughout this work, we always assume that the 



fourth-root procedure is legitimate g], and in practice it is done by setting n r = 1/4 at the end 
of the calculations. The results are 



m l} r. i 



(m x + m y ) 16*-./, 



. A\ l „2{f" , f 

(2)' 



-2a 2 ^ 2) ZRf l] ({^l})l(m 2 ) + (V ^A)+a 2 (L'( 2) +L{ 2 



j 



+ -# (4/ 3 +Pi+ 4P2) K +m d ) + -p-(-Pi- 4p 2 ) K + m y ) } , (2.3) 

J (2) J (2) 
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=/(2){ 



1 + 



16^r 2 / ( 2 2) 



3^ 



m Q B ) 



Q,B 



+ 4 ( z ( m l/) + l ( m l) + H - m I,)'>4) + « - 



+ (X o F ) + 2 £*f 11 ({^g, })/(m 2 ) ) + (V O A) 



+ a 2 (L'; 2) -L{ 2) )] +^(4l 4 -p 1 )(m u + m d ) + ^L(p l )(m x + m y )\. (2.4) 

J Z J{2) A J(2) } 



M(2) 



Here, 5(/ 2 ^ and are LO taste- violating parameters, and L'^ and L',^. are linear combinations 
of NLO taste-violating parameters. Definitions for meson masses and residue functions R and D 
can be found in Ref. [Eh. All LECs in Eqs. 



'(2) 



and (2.4) are one-loop renormalized. 
In the limit where the light valence quark masses, light sea quark masses and taste-splittings 

2 a 2 8' 

are all small compared to the strange quark mass, i.e., — ,— ZU- /W ~ £ « l we expect 

that SU(2) theory to be generated from the SU(3) theory This can be seen by expanding the 
corresponding SU(3) formulae [Q] for m 2 /(m x + m y ) and f n in powers of e. Indeed, one can check 
that the expansion has the same pattern as the SU(2) formulae. Furthermore, one can relate SU(3) 
and SU(2) LECs by comparing these two sets of formulae. First, by comparing LO meson masses 
in various taste channels, we get the relations between LO taste-violating parameters 



a 2 A^ = a 2 A#, a 1 5^ 



2 St' ( 2 ) 



2 s/ 
a o, 



Vi " °A 



,2 s:/ 



(2.5) 



Then applying Eq. (2^) in the NLO formulae of m 2 n / (m x + m y ) and f % , and comparing coefficients 
of terms m x ,m y ,mi and a 2 separately, one obtains the following relations [0]: 



/(2) 

M(2) 
Pi 
P2 

h 
I A 

L (2) 
L (2) 



1 ^ 32(2L 6 -L 4 ) 



487T 2 / 2 

1 1 km. 
-8L 8 + — -(log^+V) + 



/ 2 



167T 2 6 V te A 2 

8(2L 6 -L 4 ) + 4(2L 8 -L 5 
1 1 



— 7 (l+log^ r , 
4 , 



8L4 + 4L5 
1 



167T 2 4 

4 



;i+io g 



167T 2 4 

1 1 

16tt 2 36 

A 2 j ' 



(1+log 



L"-±A / (l+log^)--A m ,(l+log A , 



1 



1 



L'--A/( 1+log 



A 2 ^ + 2 A - (1+1 ° 8 A 2 



(2.6) 

(2.7) 

(2.8) 

(2.9) 

(2.10) 
(2.11) 

(2.12) 

(2.13) 
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where L^,L^,L^ and L% are renormalized SU(3) LECs, and L" and 11 are the NLO taste-violating 



parameters in SU(3) rSZPT. Equations ( |2.10| ) and ( |2. 1 1| ) are the same as the equations in the full 



QCD continuum case Equations ( ^8| ) and (2.9) relate the unphysical LECs in the partially 



quenched two-flavor theory to the physical LECs in the three-flavor theory. Equations ( 2.12| ) and 



( [2.13 ) give us relations between taste- violating parameters in the two-flavor and three-flavor theo- 



ries. 



3. SU(2) chiral analysis of the MILC data 

Currently, we have gauge ensembles generated from 2+1 dynamical simulations using asqtad 
staggered fermions. Lattice spacings range from 0.15fm to 0.045 fm. The light pseudoscalar mass 
and decay constant are measured with different combinations of valence and sea quark masses. 
With these data, we perform a systematic NNLO SU(2) chiral analysis by using the NNLO for- 
mulae for m 2 K and f n . Results of physical quantities are updated from the previous analysis in 



Ref. [12] 



The NNLO formulae for rrr n j{m x + m y ) and f % are obtained by combining the NLO formulae 



Eqs. ( |2.3| ) and ( |2.4j ), possible analytic NNLO terms and the continuum NNLO chiral logarithms 



provided by Bijnens and Lahde Q13TJ. The root mean square (RMS) average pion mass is used in 
NNLO chiral logarithms. In order for the SU(2) formulae to apply, we require both the valence 
and sea light quark masses to be significantly smaller than the strange quark mass. In practice, we 
used the following cutoff on our data sets: m/ < 0.2m s ,m x + m y < 0.5m s ,ma.x(m x ,m y ) < 0.3m s . 
Furthermore, in order for the continuum NNLO chiral logarithms to be applicable, we also require 
the taste-splittings between different pion states to be significantly smaller than the kaon and pion 
masses. The lattices that are at least close to satisfying these conditions are fine (a m 0.09 fm), 
superfine (a « 0.06 fm) and ultrafine (a 0.045 fm) lattices. In practice, we used superfine and 
ultrafine lattices for our central value fit, and we included fits to all three kinds of lattices to estimate 
systematic errors. One difference from the previous fit is that this time we used modified quark 
masses m — > m = m + a 2 / (2ju) in NNLO analytic terms to make the NNLO LECs scale invariant 
on the lattice, not just in the continuum. 

4. Results 

For the central fit, we used three superfine ensembles (ami,am s ) = {(0.0018,0.018), (0.0025, 
0.018), (0.0036,0.018)} and one ultrafine ensemble (am u am s ) = (0.0028,0.014). There are a 
total of 50 data points and 30 parameters with appropriate constraints. This fit has a # of 18 with 
20 degrees of freedom, giving a confidence level CL « 0.6. The volume dependence at NLO has 



been included in the fit formulae, and a small (< 0.3%) residual finite volume correction [|14J, |8J] is 
applied at the end of the calculations. 

In Fig. [I| we show the fit results for the light pseudoscalar mass and decay constant as functions 
of the sum of the quark masses (m x +m y ). The red solid curve represents the full NNLO results for 
full QCD in the continuum case, where we have set the taste-splittings and taste-violating param- 
eters to zero, extrapolated to a = linearly in a s a 2 , and set valence and sea quark masses equal. 
The continuum results through NLO and LO are shown in blue and magenta curves respectively. 

Finally, we find the physical values of the average up and down quark mass m by requir- 
ing that % has its physical mass, and then find the decay constant corresponding to this point 
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Figure 1: SU(2) chiral fits to f K (left) and ot^/(ot x + ot v ) (right). Only points with the valence quark masses 
equal (m x = m y ) are shown on the plots 



in Fig. |T(aJ With the scale r x = 0.3133(23) fm determined by HPQCD [JT5p , we obtain f n = 
130.2 ± 1.4 (+^) MeV, where the first error is statistical and the second error is systematic. This 
agrees with the PDG 2010 value f n = 130.4 ± 0.2MeV rflS]. Alternatively, one can fix the scale 



by using the SU(3) NNLO result of f n [ 17]. We then obtain 



f 2 = 123.8 ± 1.4 MeV 
h = 2.85 ±0.8l(l8;g) 
m = 3.19(4)(5)(16)MeV 



B 2 

k = 



2.91(5)(5)(14)MeV 



3.98 ±0.32 1 



+0.5 1\ 
-0.28 J 



(uu) 2 = -[281.5(3.4) (1^) (4.0)MeV] J (4.1) 
The quark masses and chiral condensate are evaluated in the MS scheme at 2 GeV. We used the 



two-loop perturbative renormalization factor [ ]18[ ] to do the conversion. Errors from perturbative 
calculations are listed as the third errors in these quantities. All the quantities agree with SU(3) 



results [17] within errors. 



5. Discussion and outlook 

In this work, we studied SU(2) rSZPT in the partially-quenched case, and we performed a 
systematic SU(2) chiral analysis for recent asqtad data in the light pseudoscalar sector. Results 
for SU(2) LECs, decay constant and chiral condensate in the chiral limit are in good agreement 
with results from an SU(3) analysis [|T^]. It can be seen that the SU(2) theory within its applicable 
region converges much faster than the SU(3) one. For the point x = 0.05 on the x-axis in Fig. [j], 
the ratio of the NNLO correction to the result through NLO is 1% both for f n and nv\j (m x + m y ). 
In contrast, the same ratio in SU(3) analysis is 3% for f n and 15% for m 2 n / (m x + m y ), respectively, 
although the large correction in the mass case is partly the result of an anomalously small NLO 
term. 
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Since the simulated strange quark masses vary slightly between different ensembles, parame- 
ters in SU(2) ZPT should also change with ensembles. In this work, we also tried to include this 
effect by using the adjustment formulae in Ref. It turns out that the fits are improved, but not 
significantly. This part may still need further investigation, hence we do not include these results 
in this work. 

In the future, a next step is to include the kaon as a heavy particle in SU(2) SZPT in order to 
study physics involving the strange quark, e.g., the kaon mass and decay constant. This method 
has been used in Ref. [|l^]. Another step is to extend the analysis to data obtained from simulations 
with HISQ fermions, where taste-violating effects are further reduced. This can be done using the 
same approach as soon as the data are available. 

We thank J. Bijnens for providing the FORTRAN code to calculate the NNLO partially- 
quenched chiral logarithms. 
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